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Abstract  (continued) 


the  statistical  behavior  of  the  signal  from  such  a  measurement  signal  has  been  developed. 
From  the  theory,  the  probability  distribution  for  the  time  intervals  between  the  zeros 
of  the  Brownian  motion  velocity  may  be  obtained. 


Continuing  research  on  this  program  will,,  focus  on  measurements  of  particles  of  known 
size  suspended  in  liquids  in  order  to  maximize  the  signal  to  noise  of  the  optical  measuring 
system,  to  refine  the  data  acquisition  methods,  and  to  verify  the  theoretical  understanding 
of  the  acquired  signal.  Following  the  demonstration  of  the  experimental  approach,  the 
technique  will  be  applied  to  measurements  of  submicron  particles  suspended  in  gas  flows. 
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1.  RESEARCH  OBJECTIVES  (STATEMENT  OF  WORK) 


The  research  objectives  for  this  program  are  outlined  in  the  Statement  of 
Work,  which  includes  the  following  tasks: 

Task  I  -  The  contractor  shall  construct  a  laser  diagnostic  apparatus  to  be 
used  for  submicron  particle  size  analysis  by  measurement  of 
Brownian  motion.  The  development  of  this  apparatus  shall  include: 

a)  Construction  and  assembly  of  relevant  optical  components; 

b)  Adaptation  of  suitable  data  acquisition  and  recording 
apparatus;  and 

c)  Development  of  appropriate  software  for  analysis  of  signals 
obtained  from  the  Brownian  motion  instrument. 

Task  II  -  Using  the  optical  apparatus  and  data  recording  and  analysis  tools 
developed  in  Task  I,  the  contractor  shall  investigate  the 
development  of  optical  heterodyne  techniques  to  be  used  in 
measurement  of  the  Brownian  motion  of  individual  particles. 
Specifically,  techniques  are  required  which  permit  measurements  of 
motion  where  the  particle  velocity  changes  in  distances  shorter 
than  one  light  wavelength. 

Task  III  -  Using  the  instrumentation  and  techniques  developed  in  Tasks  I  and 
II,  the  contractor  shall  analyze  the  statistics  of  particle  motion 
in  the  gas  to  investigate  the  particle  mass,  friction  coefficient, 
and  equilibrium  temperature.  These  studies  shall  be  conducted 
under  carefully  controlled  flow  conditions,  initially  using 
standard  size  particles.  Subsequent  development  of  the  technique 
may  utilize  a  well-characterized  flat-flame  burner. 

♦ 
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Task  IV  -  The  contractor  shall  investigate  the  measurement  of  the  refractive 
index  of  particles  in  the  Rayleigh  scattering  regime  by  analysis  of 
the  scattered  light  intensity  and  particle  mass  (size)  obtained 
using  the  Brownian  motion  sensor. 

Task  V  -  By  application  of  the  Brownian  motion  sensor  in  a  flat-flame 
burner,  the  contractor  shall  assess  the  application  of  this 
technique  for  in-situ  sizing  of  submicron  particles  in  combustion 


s  treams 


2.  RESEARCH  STATUS 


2 . 1  Laboratory  Experiments 

2.1.1  Apparatus  Development 

The  observation  of  the  Brownian  motion  of  Individual  small  particles 
poses  several  experimental  design  challenges  due  to  the  small  scattering 
cross-section  of  the  submicron  particles  combined  with  the  short  observation 
time  required  to  resolve  the  Brownian  motion.  As  a  result,  the  experiment  is 
severely  limited  by  signal  to  noise  considerations,  and  much  effort  has 
focused  on  refining  the  optical  system  to  improve  the  signal  to  noise  ratio. 

The  optical  system  which  has  been  assembled  for  this  study  is  shown 

schematically  in  Figure  2.1.  A  Spec tra-Physics  Model  165  argon  ion  laser  with 

a  single-line  power  of  approximately  2.5  watts  which  yields  a  true  TEMqq 
output  with  optimum  focusing  characteristics  is  used.  A  Spectra-Physics  Model 
589  air  spaced  etalon  is  inserted  into  the  laser  cavity  to  restrict  laser 
oscillation  to  a  single  longitudinal  mode.  This  provides  the  stable  mode 
amplitude  and  increased  coherence  length  required  in  our  optical  measurement 
system,  but  at  the  expense  of  a  75%  decrease  in  laser  power.  The  vertically 
polarized  output  of  the  laser  is  expanded  by  25x  and  collimated  to  a  beam  size 

of  approximately  3  cm.  The  beam  expansion  leads  to  tighter  focusing  and  a 

resultant  higher  power  density  at  the  sample  volume  than  can  be  achieved  with 
a  smaller  diameter  beam.  An  adjustable  iris  diaphragm  is  used  after  the  beam 
expander/collimator  to  select  the  beam  diameter. 

The  expanded  beam  is  split  into  two  beams  (b^  and  hp)  by  a  50/50  cube 
beamsplitter,  which  reflects  50%  of  the  incident  light  and  transmits  35%.  An 
increasing  intensity  difference  in  two  beams  when  focused  to  a  point  to  form 
interference  fringes  results  in  a  decreasing  clarity  of  the  fringe  pattern, 
which  is  detrimental  to  the  signal  to  noise  ratio  in  our  optical  system. 
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Figure  2.1  Schematic  of  Experimental  Optical  System  for  the  Study  o 
Brownian  Motion  of  Individual  Submicron  Particles. 


Therefore,  three  reflectors  are  used  to  decrease  the  intensity  of  the 
reflected  beam  ( 51  10%  light  loss  off  of  each  mirror)  as  well  as  to  position 

this  beam  on  a  path  parallel  to  the  transmitted  beam  (bj).  The  two  beams  of 

approximately  equal  intensity  are  then  focused  by  a  multi-element  lens  to  form 
an  interference  fringe  pattern  at  the  sample  volume.  The  sample  volume  is 
imaged  through  a  pinhole  and  then  refocused  onto  a  fast  photomultiplier 
(Type  R663)  which  views  the  light  scattered  at  90°.  An  interference  filter  is 
used  before  the  photomultiplier  to  reject  light  outside  a  1  nm  band  centered 
at  514.5  nm . 

The  measurement  volume  is  defined  by  the  intersection  of  the  laser  beam 
intensity  profile  and  the  field  of  view  provided  by  the  pinhole  of  the 

receiving  optics.  Its  geometry  is  depicted  in  Figure  2.2.  The  incident  beams 

are  at  a  7.5°  angle  with  respect  to  each  other  and  have  a  total  power  of  130 
milliwatts  with  the  iris  diaphragm  fully  open.  The  fringe  pattern  formed  from 
the  two  beams  is  easily  viewed  by  magnifying  the  measurement  volume  by  a 
factor  of  200  and  projecting  it  onto  a  screen.  The  number  of  observed  fringes 
is  consistent  with  a  calculated  4  urn  fringe  spacing  and  12  um  beam  diameter. 

A  successful  measurement  require-  the  double  coincidence  that  a  single 
particle  traverse  the  fringe  system  and  be  detected  within  the  field  of  view 
of  the  receiving  optics.  Initially,  the  measurement  volume  was  imaged  onto  a 
25  um  pinhole  with  a  magnification  of  unity,  but  low  counting  rates  were 
observed  with  such  a  small  field  of  view.  Also,  the  detection  of  scattered 
light  through  the  small  pinhole  was  highly  susceptible  to  random  mechanical 
vibrations.  A  400  jm  pinhole  was  subsequently  used  and  the  signal 
detectability  was  greatly  increased,  however,  the  measurement  volume  was  too 
large.  Extraneous  signals  from  the  interaction  of  a  single  beam  with  a 
particle  arose  because  the  field  of  view  of  400  um  was  not  restricted  to  the 
overlap  region  of  the  two  beams.  Also,  the  larger  field  of  view  increases  the 
chance  of  uetecting  multiple  particle  events.  Presently,  the  sample  volume  is 
imaged  with  a  magnification  of  =  10  onto  a  400  urn  pinhole.  This  is  equivalent 
to  a  field  of  view  of  40  um  and  renders  the  signal  less  susceptible  to 
vibra  tions . 
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Figure 


Geometry  of  the  Sample  Volume,  as  Seen  With  the 
of  View  of  25,  50,  and  WOO  ,:m  Pinholes. 


The  major  experimental  obstacle,  thus  far,  has  been  the  proper  placement 
of  the  p inhole/de  tec  tor  so  that  the  field  of  view  incorporates  the  limited 
volume  of  beam  overlap.  Several  alignment  aids  have  been  developed  to 
overcome  this  difficulty.  Spatial  overlap  of  the  two  beams  is  verified  by 
magnifying  and  viewing  the  image  of  the  focal  region,  in  a  manner  similar  to 
observing  the  fringe  pattern  described  earlier.  The  pinhole,  interference 
filter,  lens,  and  photomultiplier  are  housed  together  and  are  carefully 
baffled  to  block  out  extraneous  light.  This  unit  is  mounted  on  an  XY2 
translation  stage  to  facilitate  the  critical  detector  adjustments. 

First,  the  detector  must  be  positioned  along  the  x-axis,  i.e.,  the  axis 
along  the  scattering  direction  (See  Figure  2.1),  so  that  the  image  point  of 
the  scattered  light  is  at  the  pinhole.  This  is  best  accomplished  by  inserting 
a  mirror  at  the  focal  point  of  the  incident  beams  and  directing  the  light 
along  the  scattered  light  path.  The  detector  is  traversed  along  its  v-axis 
(perpendicular  to  the  scattered  light  path)  and,  at  the  proper  x-axis  position 
the  single  image  of  the  intense  beams  viewed  after  the  pinhole  appears  and 
disappears.  If  the  pinhole  is  not  at  the  image  point  of  the  beam  focus,  then 
traversing  the  detector  along  the  y-axis  leads  to  a  sequential  appearance  of 
radiation  from  each  beam  while  viewing  them  after  the  pinhole. 

The  optimal  y  and  2  axis  positions  of  the  detector  are  determined  by 
inserting  a  4  urn  diameter  wire  (TSI  Model  121 0— T 1.5)  precisely  into  the 
measurement  volume  focal  region  as  a  scattering  target.  The  correct  wire 
position  is  determined  by  observing  the  interference  fringes  formed  by  the 
scattered  light  and  the  primary  beams.  Then,  the  receiving  optics  are 
adjusted  for  maximum  signal. 

A  nicrocompu ter- based  data  acquisition  system  with  very  fast  time 
resolution  capability  has  been  assembled  for  this  project.  This  system  is  an 
IBM  Personal  Computer  XT  with  a  CAMAC  interface  (DSP  Technology,  Fremont, 
Calif.).  This  system  permits  considerable  flexibility  in  that  any 
CAMAC-compa tible  module  may  be  used.  Fast  time  resolution  in  data  acquisition 
is  achieved  with  the  appropriate  choice  of  CAMAC  hardware  in  combination  with 


fast  memories  for  temporary  data  storage.  The  contents  of  the  data  memories 
are  transferred  into  the  computer  memory  via  a  DMA  interface  at  data  rates  of 
approximately  400  kHz.  Large  volumes  of  data  may  then  be  stored  on  the 
system's  fixed  disk. 

Presently,  we  have  set  up  and  tested  this  system  with,  a  fast  amplifier/ 
trigger  generator  (LeCrov  Model  ol02)  and  transient  recorder  (Transiac  Model 
2001  5'  enabling  acquisition  of  up  to  32,000  (3-bit)  data  points  at  a  maximum 
rate  of  100  MHz  (10  ns  sampling  interval). 

2.1.2  Particle  Measurements 

A  necessary  step  in  the  development  of  our  interferometric  optical  system 
to  measure  individual  submicron  particles  in  a  gas  stream  is  to  verify  a 
standard  laser  velocimeter  mode  of  operation.  Various  sources  of  particles 
are  used  to  aid  in  verifying  the  "LDV”  nature  of  the  signal.  Representative 
oscilloscope  traces  of  the  scattered  light  signal  are  displayed  in  Figure  2.3. 
The  source  of  particles  in  these  traces  is  a  stream  of  H  20  droplets  directed 
through  the  sample  volume. 


Figure  2.3  Top:  Representative  Measurement  of  H -0  Particle  Traversing 
the  .Sample  Volume  (x-axis  =  s/d  i vi s  i on  ) . 

Bottom:  Single-Beam  Measurement  Taken  "nder  the  Sane 

Conditions  as  Above. 


The  upper  trace  illustrates  a  successful  measurement  of  an  individual 
particle  traversing  the  interference  fringe  pattern  formed  from  the  two 
incident  beams.  The  signal  is  nearly  100%  modulated  which  verifies  that  the 
beams  are  of  approximately  equal  Intensity  and  indicates  that  the  diameter  of 
the  measured  particle  is  less  than  or  on  the  order  of  the  calculated  4  urn 
fringe  spacing.  The  beams,  which  are  only  25%  of  their  30  mm  expanded 
diameter,  are  focused  to  a  calculated  beam  waist  of  22  um  based  on  the 
observed  number  of  fringes  in  the  trace  and  the  calculated  fringe  spacing. 

This  is  less  than  a  factor  of  two  larger  than  expected  from  diffraction 
limited  focusing  estimates  and  confirms  that  the  optics  are  of  sufficient 
quality  to  generate  the  necessary  high  power  densities  at  the  sample  volume. 

The  lower  trace  in  Figure  2.3  is  taken  with  one  beam  blocked  and, 
therefore,  represents  an  individual  particle  crossing  a  single-beam  profile, 
and  not  an  interference  fringe  pattern.  The  noise  level  in  the  single-beam 
case  is  greater  than  can  be  accounted  for  solely  by  photomultiplier  shot 
noise.  Instead,  it  is  probably  due  to  laser  noise  itself,  or  some  laser 
interaction  along  the  optical  path.  Further  characterization  and  subsequent 
suppression  of  this  noise  is  being  investigated. 

Further  verification  of  the  "LDV"  signal  has  been  achieved  by  varying  the 
velocity  of  the  droplet  stream  and  observing  the  expected  change  in  time 
between  fringe  crossings  by  the  particle  (At  ■  X/2nVsin  9/2,  where  V  is  the 
velocity  component  in  the  90°  scattering  direction).  Also,  successful 
measurements  have  been  recorded  for  individual  particles  of  Ti02  in  a  gas 
stream  formed  from  the  chemical  reaction: 


*  2H2°<g) 


TiO_  /  \  +  4HC1,  \ 
2(s)  (g) 


Although  the  size  distribution  of  the  Ti02  particles  generated  by  this 
reaction  has  not  been  well  characterized,  similar  formation  experiments  have 
yielded  Ti02  size  distributions  between  0.5  and  0.05  um  in  diameter.  Future 
plans  include  direct  characterization  of  the  Ti02  particles  formed  in  the 
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laboratory  and  preparation  of  other  sources  of  submicron  particles  in  a  gas 
stream. 

A  complementary  and  perhaps  more  exact  means  of  calibrating  and 
characterizing  the  optical  system  involves  the  measurement  of  particles 
suspended  in  liquid  solutions.  To  this  end,  a  particle  flow  cell  has  been 
designed  and  built  in  the  laboratory  which  enables  us  to  measure  the  ”LDV” 
signal  from  particles  of  known  diameter  with  the  added  advantage  of  well 
characterized  flow  velocities.  These  calibration  particles  are  commercially 
available  (Coulter  Electronics)  in  a  variety  of  sizes  in  the  submicron  range. 
Preliminary  measurements  of  0.261  urn  diameter  latex  particles  have  been 
recorded,  however,  the  added  alignment  difficulties  associated  with  the  flow 
cell  requires  further  refinement  of  the  alignment  procedures.  These 
refinements  are  currently  in  progress. 

2.2  Behavior  of  the  Brownian  Motion  Signal 
2.2.1  Introduction 

In  the  last  report,1  a  Monte  Carlo  simulation  of  the  signal  from  the 
Brownian  motion  sensor  was  carried  out.  The  result  of  interest,  namely,  the 
time  between  zero  crossings  of  the  Brownian  velocity,  was  found  to  be 
dependent  on  the  step  size  At  of  the  simulation.  This  dependence  is  actually 
not  surprising  because  At  represents  essentially  the  bandwidth  of  the  actual 
signal  detection  system,  which,  being  a  practical  device,  should  have  a  finite 
bandwidth.  Since  the  Brownian  motion  is  a  statistical  process,  the  behavior 
of  any  signal  derived  from  the  Brownian  motion  is  statistical,  and  the 
spectral  properties  of  the  signal  would  depend  on  the  bandwidth  of  the  signal 
detection  system.  These  ideas  are  further  explored  here. 

In  the  following,  the  characteristics  of  a  Brownian  motion  detection 
system  are  first  described;  the  most  important  result  being  that  the  only 
meaningful  information  carried  by  the  signal  is  contained  in  the  time  between 
the  extrema,  r,  of  the  signal.  Then  the  theory  on  the  statistical 
distribution  of  the  values  of  t  is  presented. 
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The  configuration  of  a  Laser  Doppler  Brownian  motion  detection  system  has 
been  discussed  in  Reference  1  and  elsewhere  in  this  report.  For  a  perfect 
detection  system  as  a  whole,  the  transfer  function  is  shown  in  Figure  2. A.  By 
transfer  function,  we  mean  the  relationship  between  the  input  quantity  of 
Interest,  namely  the  position  of  the  particle  in  Brownian  motion,  and  the 
output  signal  which  is  the  beat  amplitude  of  the  two  fringe  forming  laser 
beams  as  seen  by  a  square- law  detector. 


ik  •£  ik^  *£ 

signal  '  <|Re(^e  +  EL,e  )|>  (D 

'1  +  cos 

In  the  normal  LDV  application,  the  particle  traverses  the  illuminated  volume 
with  a  constant  velocity  (from  A  to  B  in  Figure  2. A).  If  u  is  the  velocity 
component  in  the  direction  of  k j -k ? , 

AJc*_r  3  ut,  (2) 

the  signal  will  be  modulated  fully  (A'B*)  as 

signal  -1  +  cos  (Akut)  (3) 

and  the  velocity  may  be  obtained  from  the  period  t  of  the  modulation, 

u  ■  2ir/Ak  x  (A) 

Note  that  implicit  in  this  velocity  measurement  is  a  length  scale  (Ak)"1  so 
that  a  time  measurement  (t)  can  yield  a  velocity. 

In  the  detection  of  Brownian  motion,  however,  u  is  a  statistical  quantity 
which  changes  rapidly  in  time.  In  particular,  the  particle  may  change 


/  s 


velocity  appreciably  (or  may  even  reverse  its  direction)  before  it  traverses  a 
"fringe”  of  the  transfer  function.  Therefore  the  signal  is  no  longer  fully 
modulated  and  (Ak)"1  ceases  to  be  a  relevant  length  scale.  Because  of  the 
nonlinear  nature  of  the  transfer  function,  the  amplitude  of  the  signal  will 
not  carry  meaningful  information,  e.g.,  when  the  particle  oscillates  about  CD 
in  Figure  2.4,  the  signal  response  C'D*  is  quite  different  from  that  (E'F')  of 
a  similar  oscillation  EF. 

The  only  meaningful  information  contained  in  the  signal  is  the  time 
between  the  extrema  of  the  signal  modulation.  Each  extremum  of  the  signal 
represents  a  reversal  in  direction  of  the  particle,  and  therefore,  a  zero  of 
the  Brownian  velocity  u.  In  order  to  interpret  the  signal,  it  is  necessary  to 
examine  the  statistical  distribution  of  the  zeros  of  u. 


2.2.3  The  S tatis tical  Nature  of  the  Brownian  Velocil 


The  Brownian  velocity  of  a  particle  in  thermal  equilibrium  in  a 
background  fluid  is  described  by  the  Langevln  equation.  For  the  observed 
velocity  component  u  ( in  the  direction  of  Ak) ,  the  Langevin  equation  is 


Su 


f 


(5) 


where  8_1  is  a  relaxation  time  governed  by  the  mass  m  and  the  mobility  t  of 
the  particle  in  the  fluid. 
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The  statistical  nature  of  u  is  due  to  the  random  driving  force  f,  which  is  the 
result  of  the  impulse  arising  from  each  collision  of  the  particle  and  a 
molecule  of  the  background  fluid.  Because  these  collisions  occur  randomly,  f 
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may  be  considered  to  be  a  white  noise  function  with  a  constant  spectral 
density  Sf, 


<  f(t)f(t*r)  >  -  Sf$(t) 


(7) 
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The  time  scale  associated  with  the  fluctuation  of  f  is  of  the  order  of 
the  time  between  the  molecular  collisions  tc,  which,  for  a  100  A  particle  in 
a  gas  medium  at  STP  is  approximately  10”11  s.  The  time  response  of  practical 
instruments  used  to  observe  the  Brownian  motion  velocity  is  usually  much 
slower  than  this  collision  time.  The  observations,  are,  therefore,  the  result 
of  averaging  u  over  a  response  time  scale  At  associated  with  the  instrument, 
with  At  >  tc.  The  observed  value  of  u  is  then  a  random  variable,  and,  since 
i t  is  the  aggregate  of  many  similar  events,  it  would  have  a  Gaussian 
distribution  by  the  central  limit  theorem.  If  the  particle  has  a  velocity 
uQ  at  one  instant  of  time,  the  conditional  probability  distribution  of  u  at 
a  time  At  later,  (so  that  the  value  of  u  is  the  aggregate  of  the  sum  of  all 
the  collisions  in  At),  is  given  by. 


p  (u,uq ;A t) 


exp  [-(u-u  h)  /a 
a/TT  r  1  o 


2kT(l-h  ) 


(8) 


(9) 


Note  that  the  distribution  Is  parameterized  by  the  "memory"  factor  h 


h  =  e 


-8  At 


(10) 
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and  the  width  a.  The  observed  velocity  u( t)  may  be  viewed  as  observations  at 
discrete  points  in  time,  tn  *  nAt.  When  the  time  step  At  Is  large  compared 
to  the  relaxation  time  (8At  >  1),  the  velocities  u^  are  independent  of  each 
other,  and  each  obeys  a  Maxwellian  distribution  at  temperature  T.  In  the 
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Brownian  motion  sensor  described  here,  since  the  velocity  amplitude  is  not 
measured,  all  the  information  is  contained  in  the  measurement  of  the  value  of 
8.  Therefore,  the  instrument  response  time  At  should  be  such  that  8At  <  1. 
After  8  is  determined,  the  particle  physical  properties  (c,m)  may  be 
determined  from  Eq.  (6). 

The  Distribution  of  the  Zeros  of  the  Velocity 

The  Brownian  motion  sensor  measures  the  time  between  successive  zeros  of 
the  Brownian  velocity  u  (in  the  direction  of  Ak) .  The  time  history  of  the 
velocity  values,  u( t) ,  as  observed  by  an  instrument  with  an  integration  time 
At,  may  be  thought  of  as  a  time  series  un  at  time  nAt.  From  Eq.  (8)  the 
values  of  each  Uj,  are  governed  by  the  probability  distribution 

P  <v  Vi*  40  ■  a7«  exp  t  "  Vih)2/,21  <n) 

where  the  un_^  appears  as  a  parameter.  We  want  to  determine  the  probability 
P(N)  of  obtaining  a  time  between  two  successive  zeros  of  NAt. 

The  problem  of  determining  the  distribution  of  zeros  of  a  random  variable 
has  been  studied  extensively  in  the  context  of  shot  noise  in  electronic 
devices,  and  in  queueing  theory.  The  classic  review  of  the  subject  is  given 
in  Reference  3.  We  shall  derive  the  expression  for  P(N)  independently  here. 

To  connect  our  result  to  that  in  the  literature  of  noise  current,  the  random 
driving  force  f  in  Eq.  (5)  may  be  considered  as  the  randomly  arriving 
electron,  and  the  observed  "current"  u  is  the  result  of  two  transfer  functions 
(or  filters)  in  series.  The  first  transfer  function  is  that  associated  with 
the  Langevin  equation.  Formally 
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where  f  is  the  Fourier  Transform  of  f.  Physically,  this  transfer  function  is 
due  to  the  averaging  of  the  effect  of  f  by  the  damping  of  the  particle  in  the 
fluid.  The  second  transfer  function  is  due  to  the  averaging  by  the 
instrument.  Therefore  the  overall  transfer  function  is 


u(t>  =  ^/  df  7 


f  e  doj 
( S  tioi^  2ir 


It  should  be  noted  that  no  closed  form  solutions  for  P(N)  have  been 
worked  out,  except  for  rather  trivial  cases.  The  difficulty  was  attributed3 
to  the  lack  of  an  analytical  solution  to  the  N  dimensional  integral  involved. 

We  shall  derive  the  expression  for  P(N)  here.  The  process  of  obtaining  a 
successive  2ero  in  N  steps  may  be  described  as  the  following.  At  time  zero 
(n  =  0),  the  particle  has  a  very  small  velocity  6u,  which,  within  the 
amplitude  resolution  of  the  instrument,  may  be  considered  as  zero.  In  the 
following  argument,  6u  is  taken  as  positive.  The  argument  will  be  exactly  the 
same  if  5u  is  negative.  We  shall  define  a  zero-crossing  event  at  t  =  NAt, 

(n  =  N),  as  an  event  with  velocity  greater  than  zero  at  t  =  (N  -  1)  At, 

(n  =  N  -  1),  and  velocity  less  than  or  equal  to  zero  at  t  =  NAt.  The  time 
resolution  of  this  definition  of  zero-crossing  is  At,  and  so  the  meaning  of 
such  an  event  for  N  =  1  is  not  precise.  Then  the  condition  for  occurrence  of 
a  successive  zero  in  N  time  steps  is  to  have  the  velocities  be  positive  for  n 
*  1,  ...,  N-l;  and  the  velocity  be  negative  or  equal  to  zero  for  n=N. 

Therefore 


P(N )  =  /  dUj^  p  (u^,  0;  At)  /  du2  p  (u2 ,  ^  ;  At) 


J  dUn-l  P  (Vl*  Un-2  ’  At)  /  dUn  P  (V  Vi 5  At) 

o 


2-14 


Using  the  formula  for  p  in  Eq.  (11),  and  changing  the  variable  from  u  to 
£  =  u/o,  the  above  probability  becomes 


-r2  . 
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Note  that  a  does  not  appear  in  this  expression.  This  is  because  the 
zero-crossing  time  pertains  to  the  zeros  of  the  velocity,  and  is  therefore 
independent  of  the  velocity  scale.  (The  value  o  is  a  velocity  scale  measuring 
the  width  of  the  velocity  distribution). 

If  the  integration  time  At  is  large  (3At  »  1),  h  will  approach  zero  and 
Eq.  (15)  reduces  to  a  binomial  process 


P  (N )  =  (i) 


The  average  value  of  N  is 


N  =  I  NP(N)  -  •  \  -  -2- 

N=1  N  =  1  2N 


|x=l  N-l 


In  this  case,  the  data  values  corresponding  to  the  successive  zeroes  of  the 
velocity  become  a  completely  random  sequence,  and  do  not  carry  information 
about  h. 

For  small  values  of  h,  the  multiple  integral  may  be  worked  out  as 
follows:  An  operator  In  is  defined  as 
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The  quantity  £n_i  enters  as  a  parameter  In  the  above  definition.  Then 


I  (1)  =  \  (1  +  r  C  ,h)  +  0(h2) 

n  2  »7r  n— l 
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The  expression  for  P(N),  Eq.  (15),  after  evaluating  the  last  integral  on  the 
right  hand  side,  becomes 


P(N>  -  \  Ij  I2  •••  IN_X  U  +  /§  +  0(h2) 


The  expression  (19)  and  (20)  may  be  used  repeatedly  In  (21).  The  first  few 
expressiono  are: 


P(N)  -  [j)  I ^  ...  IN_2  [U  -  ^  )  +  ^N_2n]  +  °^h  ^ 
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P(N)  =  (yj  [l  +  (N-3)  ^  ]  +  0(h);  N  >  2 


The  above  expression  is  valid  for  N  _>  2.  For  N  =  1 ,  since  the  velocity  at 
n  =  0  is  precondi tioned  to  be  =0, 
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Compared  to  the  binomial  distribution  (P(N)  =  2-N ) ,  which  is  the  limit 
for  h  =  0,  P(N)  decreases  with  h  for  N  <  3  and  increases  with  h  for  N  >  3. 

This  result  is  illustrated  in  Figure  2.5.  Since  the  velocity  at  n  =  0  is 
zero,  the  velocity  uj  at  n  =  1  is  sampled  from  a  Gaussian  centered  at  zero. 

The  value  of  P(N  =  1)  is  the  probability  that  this  velocity  is  negative,  and 
P(l)  is  therefore  1/2.  Say  the  velocity  at  n  =  1  is  u2(>  0),  then  the 
velocity  u2  at  n  =  2  is  sampled  from  a  Gaussian  centered  at  u^.  The 

probability  of  u2  being  negative  is  the  shaded  area  in  Figure  2  under  the 

Gaussian  at  n  =  2,  and  this  value  is  <  1/2  for  h  >  0.  Therefore  P(2) 

<  (1/2)  .  For  N  =  3,  the  probability  of  obtaining  a  consecutive  zero  here  is 

represented  by  the  shaded  area  under  the  Gaussian  centered  at  u2h  at  n  =  3 
(Figure  2.5).  This  area  has  to  be  summed  over  all  the  realizations  of  the 
possible  values  of  positive  u2,  which  is  sampled  from  the  Gaussian  'entered  at 
u2h,  which  again  has  to  be  summed  over  all  the  realizations  of  uj.  Although 

the  shaded  area  at  n  =  3  is  less  than  1/2,  there  are  "more"  realizations  than 

the  h  =  0  case  because  u2  is  sampled  from  the  shifted  Gaussian.  These  two 
effects  balance  each  other  to  0(h2),  and  P(3)  is  equal  to  (1/2) 3  +  0(h2). 

For  N  >  3,  P(N)  is  given  by  Eq.  (22),  and  is  greater  than  (1/2)N.  The 
values  of  P(N)  are  plotted  versus  N  in  Figures  2.6  and  2.7. 

I t  is  not  possible  to  obtain  a  closed  form  solution  of  P(N)  for  arbitrary 

h  (h  <  0  <  1  corresponding  to  0  <  3At  <  °°)  because  of  the  difficulty  in 

evaluating  the  multiple  integral  in  Eq .  (11).  Direct  numerical  calculation  of 
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7  Probability  of  Obtaining  a  Consecutive  Zero  of  the 
Brownian  Velocity  in  N  Steps  (Semi-Log  Scale). 
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the  Integral  is  not  feasible  for  any  large  N,  say  N  >  5.  Monte  Carlo 
evaluation  of  the  Integral  would  be  equivalent  to  our  earlier  simulation  of 
the  Brownian  motion.1  The  results  from  Reference  1  will  be  reexamined  in  this 
context  here. 

The  general  behavior  of  the  mean  number  of  steps  N  between  consecutive 
zero-crossings  of  the  Brownian  velocity  is  shown  in  Figure  2.8.  For  h  ■  0,  N 
is  2  (from  Eq.  (17)).  For  small  h, 

00  N 

N  -  l  N  (■=•)  ( 1  +  (N-3)  •££)  +  0(h2) 

N-l 

-i  <i  -  ?>  +  i  (4)  (?) +  o^2’ 

N-l  2W  N-l  2N 

N  -  0(h2)  (24) 

In  the  last  expression,  we  made  use  of  £(N/2^)  -  2  and  £(N2/2n)  -  6. 

Therefore  N  *  0(h2)  for  small  h.  For  large  relaxation  time  (8  -*•  0,  8At  0, 
and  h  +  1),  the  particle,  once  it  acquires  a  positive  (or  negative)  velocity, 
will  retain  that  velocity  for  a  large  time.  Consequently  N  +  ®. 

The  results  of  the  Monte  Carlo  simulation  in  Reference  1  have  been 
replotted  in  the  context  of  the  present  analysis  in  Figure  2.9,  and  confirm 
the  theory.  The  results  in  Figures  2.8  and  2.9  show  that  to  obtain  h,  (and 
therefore  the  values  of  8,  which  would  lead  to  the  values  for  €/n>)>  it  is 
necessary  to  have  a  time  resolution  of  At  such  that  8At  <  1  (or  h  ■  e”^At 
close  to  1 ) . 

2.2.4  Conclusion 

A  theoretical  analysis  of  the  probability  of  obtaining  a  time  t  between 
consecutive  zero-crossings  of  the  Brownian  velocity  u  has  been  presented. 
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Figure  2.8  General  Behavior  of  the  Expectation  Value  of  the  Number 
of  Time  Steps  Between  Consecutive  Zero-Crossings  of  the 
Brownian  Velocity,  as  a  Function  of  the  Memory  Parameter 
h  ■  exp  (-8At). 


h  ■  exp  (-6At) 

Figure  2.9  Monte  Carlo  Simulation  Result  of  the  Mean  Time  Step  (N) 

Between  Zero-Crossings  of  the  Velocity,  (a)  as  a  Function 
of  BAt;  (b)  as  a  Function  of  h  -  e“8At. 


Since  u  Is  a  random  variable,  the  probability  P(t)  is  necessarily  a  function 
of  the  integration  time  At  of  the  instrument  used  to  observe  u.  The  analysis 
has  been  carried  out  in  terms  of  discretizing  the  time  in  steps  of  At, 
t„  -  nAt,  and  the  resulting  probability  P(N),  (N  -  T/At),  has  been 
derived.  The  values  of  P(N)  are  found  to  be  dependent  on  the  "memory”  h  (h  h 
exp  (-BAt))  of  the  Brownian  velocity,  where  B”1  is  the  relaxation  time.  For 
very  short  memory  (h  <  1),  the  probability  P(N)  tends  to  the  binomial  limit 
P(N)  ■  1/2**.  The  expectation  value  N  (  =  £  NP(N)  assumes  the  value  of  2  in 

the  binomial  limit  and  2  +  0(h2)  for  small  h.  For  large  relaxation  time  (6  ♦ 
0),  N  -*•  ».  To  obtain  accurate  information  from  the  measurement  of  N  about  h, 
and  therefore,  the  value  of  S  which  is  related  to  the  physical  properties  of 
the  particle  suspended  in  the  fluid,  it  is  necessary  to  have  BAt  <  1.  Since 
the  result  is  a  function  of  the  integration  time  At,  which  is  not  a  precisely 
defined  quantity  for  the  measuring  instrument,  the  Brownian  motion  sensor  must 
be  a  calibrated  instrument. 

2.3  Research  Plan 

The  research  plan  for  the  third  year  of  the  program  will  continue  to 
emphasize  experimental  validation  of  the  measurement  concept.  Specifically, 
we  anticipate  performing  the  following  tasks: 

Task  1. 

We  will  continue  the  measurements  of  particles  of  known  size  suspended  in 
liquids  in  order  to  maximize  the  signal  to  noise  of  our  optical  measuring 
system,  to  refine  our  data  acquisition  methods,  and  to  verify  our  theoretical 
understanding  of  the  acquired  signal.  Data  from  these  commercially  available 
particles  in  the  0.046  ym  to  0.5  ym  size  range  should  help  establish  the 
limitations  of  our  apparatus.  Analysis  of  the  signal  from  the  different  size 
particles,  when  compared  to  the  Monte  Carlo  simulations,  should  establish 
whether  the  signal  fluctuations  exhibit  the  expected  behavior  for  particles  of 
different  size,  and  should  provide  a  crucial  test  for  the  measurement  theory. 
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Task  2. 
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If  Che  tests  with  liquid  suspensions  conclude  with  favorable  results,  we 
will  proceed  with  measurements  of  submicron  particles  suspended  in  gas  flows. 
These  measurements  will  proceed  in  two  steps: 

a.  The  well  defined  liquid-suspended  particles  described  in  Task  1  will 
be  dispersed  into  a  gas  stream.  The  results  from  such  monodlsperse 
particulates  will  be  compared  to  the  measurements  obtained  in  liquid 
suspensions. 

b.  We  are  currently  researching  methods  of  chemical  generation  of  small 
particulates  from  the  gas  phase.  Likely  methods  Include  the  Ti02 
generation  described  in  this  report  and  photochemical  particulation 
induced  by  a  pulsed  N2  laser  in  a  CS2  vapor.  The  size  distribution 
results  obtained  from  these  gas-borne  particulates  in  our  Laser 
Doppler  measurement  system  will  be  compared  to  distributions 
obtained  via  an  independent  method  such  as  electron  microscopy.  The 
measurements  of  submicron  particles  suspended  in  well  controlled  gas 
flows  should  be  critical  in  establishing  the  capabilities  and 
limitations  of  the  technique  as  a  particle  diagnostic. 
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3.  PUBLICATIONS 

No  results  to  date  from  this  program  have  been  published  in  technical 
journals.  We  anticipate  two  publications  in  the  coming  year,  one  reporting 
the  theory  and  experimental  validation  of  the  Brownian  motion  measurement 
approach  (to  be  submitted  to  Journal  of  Applied  Physics  or  Physics  of  Fluids) 
and  a  second  paper  describing  a  measurement  application  ( to  be  submitted  to 
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6.  NEW  DISCOVERIES  OR  INVENTIONS 


The  experimental  approach  under  study  in  this  program,  i.e., 
determination  of  particle  size  by  measurement  of  the  Brownian  motion  of 
individual  particles,  is  a  new  approach  which  has  not  previously  been 
demonstrated.  We  anticipate  that  this  approach,  if  successful,  will  permit 
the  development  of  research  instrumentation  for  the  study  of  the  growth  of 
submicron  particles,  based  on  new  concepts  in  measurement  and  data  analysis. 
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